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ABSTRACT 

Quantum gravity is investigated in the limit of a large number of space-time dimensions, using 
as an ultraviolet regularization the simplicial lattice path integral formulation. In the weak field 
limit the appropriate expansion parameter is determined to be 1/d. For the case of a simplicial 
lattice dual to a hypercube, the critical point is found at kc/X = 1/d (with k = I/SttG) separating 
a weak coupling from a strong coupling phase, and with 2d? degenerate zero modes at kc- The 
strong coupling, large G, phase is then investigated by analyzing the general structure of the 
strong coupling expansion in the large d limit. Dominant contributions to the curvature correlation 
functions are described by large closed random polygonal surfaces, for which excluded volume effects 
can be neglected at large d, and whose geometry we argue can be approximated by unconstrained 
random surfaces in this limit. In large dimensions the gravitational correlation length is then found 
to behave as | log(A;c — k)\^?^, implying for the universal gravitational critical exponent the value 
= 0 at d = oo. 
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1 Introduction 


The lack of perturbative renormalizability for quantum gravitation in physical dimensions 
[1, 2, 3, 4, 5, 6] has brought to the forefront the need to develop field theoretic approximation 
schemes that do not rely on the assumption of weak gravitational fields, and which are sophisti¬ 
cated enough to deal with the rich physical structure of non-renormalizable theories [7, 8, 9, 10, 11]. 
The hope is that more powerful covariant methods, better suited to the non-perturbative regime, 
will eventually shed some light on the elusive long distance properties of quantum gravitation, 
which could ultimately have a bearing on a number of long standing and fundamental issues, such 
as the short distance nature of space-time, the emergence of the semiclassical limit and the problem 
of large-scale quantum cosmology [12, 13, 14, 15]. 

Approaches based on the simplicial lattice formulation for gravity [16, 17, 18, 19, 20, 21, 22, 
23, 24], the 2 -|- e expansion [8, 25, 26, 27, 28] and approximate renormalization group methods 
based on Wilson’s momentum slicing technique [29, 30] have suggested the existence of a nontrivial 
ultraviolet fixed point in and around four dimensions, separating a weakly coupled (but physically 
un-attractive) phase from a strongly coupled one, the latter phase being characterized by a finite 
invariant correlation length, and close to smooth geometries at large distances. Substantial un¬ 
certainties remain in each of the three approximation methods mentioned above, both about the 
results themselves and their relationship to each other, but also regarding their ultimate physical 
significance and how they might relate to physical gravitational phenomena, and both early and 
late time cosmology. It would be clearly desirable if one could find a limit in the quantum gravity 
case where non-perturbative aspects of the theory could be fully explored by covariant analytical 
means. In the significantly simpler Yang-Mills case the evidence so far is that the lattice is the 
only reliable non-perturbative method, capable of producing reasonably unambiguous quantitative 
results, within a controlled approximation based on the zero lattice spacing limit. It will be this 
method that will be therefore the focus of our work. 

In this paper we study a set of approximation methods based on an expansion in the inverse 
number of dimensions. Increasing the number of space-time dimensions above four only worsens the 
renormalizability problem, which implies that the need for a non-perturbative approach, such as 
the lattice one, becomes even more acute. The so-called 1 jd expansion was originally developed for 
statistical mechanics systems, and later extended to the study of quantum field theory, where it has 
since met with a number of considerable successes, including an understanding of triviality for scalar 
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field theories above four dimensions (which, incidentally, are not perturbatively renormalizable for 
any d > 4). The above expansion is known to be intimately tied up with the mean held theory 
treatment of quantum mechanical systems, but not necessarily equivalent to it (as was already noted 
in the gauge theory case), and exploits the fact that in large dimensions each point is typically 
surrounded by many neighbors, whose action can then be either treated exactly, or included as 
some sort of local average. For classical spin systems at hnite temperature, the 1/d expansion was 
originally developed in [31, 32, 33] by examining the structure of the high temperature expansion. 

In many ways the 1/d expansion is similar to the very successful 1/iV expansion for statistical 
mechanics systems (the 0{N) vector model being one thoroughly explored and well-understood 
example [34, 35]) and SU{N) gauge theories, where it leads to the planar diagram approximation 
[36, 37] and the many phenomenological successes that follow from it. In the gravitational case it is 
less obvious how to attach color degrees (or any other internal degree) of freedom to the graviton, 
so this particular avenue seems unfruitful at the moment. 

In this paper we will study large-dimensional pure gravitation, without any matter fields, which 
could then be added at a later stage. We recall here that for pure gravity in d dimensions there are 
d(d-|-l)/2 independent components of the metric, and the same number of algebraically independent 
components of the Ricci tensor. The contracted Bianchi identities then reduce the count by d, and 
so does general coordinate invariance, leaving d(d-|- l)/2 — d — d = d(d — 3)/2 physical gravitational 
degrees of freedom in d dimensions. As a result, the number of physical degrees of freedom of the 
gravitational field grows rather rapidly (quadratically) with the number of dimensions. 

The paper is organized as follows. In Section 2 we discuss the machinery of the 1/d expansion 
for the lattice theory of gravity based on Regge’s simplicial construction. The action simplifies 
considerably in the large d limit, and we are able to exhibit the location of the critical point in the 
variable k = I/SttG, at least in the weak field limit, as well as the nature of the excitation spectrum 
around it. In Section 3 we follow a complementary route to the large d limit, where we perform a 
simultaneous 1/d and strong coupling (small k) expansion. Since the strong coupling expansion for 
simplicial lattice gravity has not been discussed before in the literature, we will present here some 
general aspects of it. We then show how the relevant (in the long distance limit) critical behavior 
can be extracted from the strong coupling expansion by analyzing the geometric structure of its 
dominant terms. In Section 4 we provide some contact with results obtained in the continuum in 
and above d = 4, and compare and contrast with what has been found in the previous two sections 
from the simplicial lattice theory. Appendix A contains a brief summary of the large d limit for 
scalar lattice field theories, while Appendix B discusses some results relevant to non-Abelian gauge 
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fields on the lattice. 


2 Expansion in Inverse Powers of the Dimension 


Our first concern will be an approximate evaluation, in the large d limit, of an appropriately 
discretized form of the continuum Euclidean functional integral for pure gravity without matter, 
which we write here as 

Zcont = J [dgf,u] exp ^-A J d'^Xy/g + J d'^Xy/gR^ . (2.1) 

In the following we will therefore first address the key issue of precisely what type of terms in the 
discrete action, based on the simplicial lattice formulation [16], become dominant in this limit. 

2.1 General formulae in d dimensions 


We will consider here a general simplicial lattice in d dimensions, made out of a collection of flat 
d-simplices glued together at their common faces so as to constitute a triangulation of a smooth 
continuum manifold, such as the d-torus or the surface of a sphere. Each simplex is endowed with 
d + 1 vertices, and its geometry is completely specified by assigning the lengths of its d{d + l)/2 
edges. We will label the vertices by 1, 2, 3,... d + 1 and denote the square edge lengths by I 12 = ^21 > 
• •• l-id+i- vertices of the simplex can be specified by a set of vectors vi = 0, V 2 , ... v^+i. The 
matrix 

gij = < Vi+i|vj+i > , (2.2) 


with 1 < i, j < d, is positive definite, and, in terms of the edge lengths Uj = |vj — Vj|, it is given 
by 

dij — 2 (^M+1 • (2-3) 

The volume of a d-simplex is then given by the d-dimensional generalization of the well-known 
formula for a tetrahedron 

Vd = ^y^det^ . (2.4) 

An equivalent form can be given in terms of a determinant of a (d -|- 2) x (d -|- 2) matrix, 

,1/2 


Vd = 


(- 1 ) ^ 
d! 2‘^/2 


0 

1 

1 

1 

0 

/2 

tl2 

1 

^21 

0 

1 

llx 

/2 

‘32 

1 

/2 

V+1,1 

/2 

V+1,2 • • • 


(2.5) 
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Then the dihedral angle in a d-dimensional simplex of volume Vd, between faces of volume Vd-i 


and is obtained from 


sinOd = 


d VdVd-2 


d -1 Vd-i yj.i 

In the equilateral case we record here the particularly simple result for the volume of a simplex 


Vd = 


1 / d + 1 
2'^ 


and for the dihedral angle 


cos 6d = 


The d-dimensional Euclidean lattice action, involving cosmological constant and scalar curvature 


terms, is then given by 


/(/ 2 ) = xJ2Vd - kY, SdVd-2 , 


and appears in the partition function as 


Z{X,k) = J[df] 


( 2 . 10 ) 


2.2 Weak field expansion 

The above formulae for volumes and angles are quite complicated in the general case, and therefore 
of limited use in large dimensions. The next step consists in expanding them out in terms of small 


edge length variations. 


fi. = 0^ + Sfi- 


( 2 . 11 ) 


We will set for convenience from now on 5 ifj = etj. Unless stated otherwise, we will be considering 
the expansion about the equilateral case, and set = 1 (we will later relax this last restriction). 
Furthermore one has the well known expansion for determinants 

det(l + M) = 

= l+trM+— (trM)^ - trM^ +— {tr Mf - 3 tr M tr + 2 tr + ... 

( 2 . 12 ) 


One can then re-write the expression in Eq. (2.5) for the volume of a d-simplex as 

{ ll^^ 

= -nw- ' 

and expanding out to quadratic order one finds 

a/— det M 2 = a/s -|- —j= ei 2 + ... + €12 ei3 + • • • — 0 + ■ ■ ■ 


(2.13) 


(2.14) 
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+ +J^ei2ei3 + ... -^€12634 + .+ (2-15) 


and for general d 

, . -, %/=*= det Md = 1 + , ei2 + 

Vd +1 d+l 


(dTlF ''' " JdTTr " 2(d+l)2 + • • • + o(e')- 


(2.16) 


For large d the last expression simplifies to 


^ 'J ± det Md — 1 + ^ (ci 2 ^ (^12 ^13 + • • •) ~ ^ (^12 £34 + • • •) ~ 2 ^^12 + •••) + 

(2.17) 

Here the terms ei 2 €43 refer to two edges sharing a common vertex, whereas the terms 642 €34 denote 
terms with opposite edges, not sharing a common vertex. 

As a result, the volume term appearing in the d-dimensional Euclidean lattice action of Eq. (2.9), 


becomes 


d^oo d!2^/2 


{1 + ^(^12 + ...) + ^(642 643 + ...)- ^(612634 + ...)- 2 (^12 + ' ' O + •••|) 

(2.18) 


or, equivalently, ordering the terms in powers of 1/d, 


d\ 2'^/2 


1 ~ X ^42 + • • • + ”7(^12 + • • • + ^42 643 + . . .) + 


(2.19) 


To leading order, it involves a lattice sum over all squared edge length deviations. Note that the 
terms linear in e (the so called tadpole terms in the continuum), which would have required a shift 
in the ground state value of e for a non-vanishing cosmological constant A, vanish to leading order 
in 1/d. The full volume term X^Vd appearing in the action can then be easily written down using 
the above expressions. 


Next one needs to expand the dihedral angle. In the equilateral case one has for the dihedral 


angle 


9d = arcsin ■ 


d d- 


TT 1 


d Hd? 


+ ... , 


( 2 . 20 ) 


which will require four simplices to meet on a hinge, to give a deficit angle of 27r — 4 x ^ 0 in large 

dimensions. One notes that in large dimensions the simplices look locally (i.e. at a vertex) more 
like hypercubes. Several d-dimensional simplices will meet on a (d—2)-dimensional hinge, sharing a 
common face of dimension d—1 between adjacent simplices. Each simplex has (d—2)(d —1)/2 edges 
“on” the hinge, some more edges are then situated on the two “interfaces” between neighboring 


simplices meeting at the hinge, and finally one edge lies “opposite” to the hinge in question. In 
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two dimensions one finds for the dihedral angle at vertex 1, to quadratic order, 


whereas in three dimensions one has for the dihedral angle at edge 12, to the same order, 

^ . 2^/2 1 1 , ,1 
03 = arcsm — -h ei 2 - (eis + cm + 623 + ^24) + ;^ £34 

+ + '■* + 4) - 5^‘12(£1S + £14 + £23 + £24) 

— -— 7 =(ei 3 e 24 + ^14623) + -— 7 =(ei 3 ei 4 + £23^24) + /- (^13 €23 + ^14624) 

AyZ AyZ SbyZ 


+ 41 £12 £34 - 4?! 4s + £14 + £23 + £24) £34 - £34 


( 2 . 22 ) 


In the general d-dimensional case the expansion coefficients for the dihedral angle at the hinge 
labeled by 1, 2,... d — 1 are given by the following expressions (as well as their large d limit) 


dVdP^ 

d- 1 

d- 1 

/ 1 ^d,d+l 


2{d^ -2d^ -d+l) 2 

d2(d2 - 1)3/2 ^12 

(d^ — 2d — 2)(d — 1)^ 2 
2d2(d2 - 1)3/2 

id-If .2 

2(d2 - 1)3/2 ^d,d+l 

2{d? - 4d2 - d + 2) 

d2(d2 - 1)3/2 ^12 £13 

4(2d2 - 1) 

“d2(d2 - 1)3/2 ^12 £34 

(d-l)(d3-2d2 + d + 2) 

d2(d2 - 1)3/2 ^12 £l,d+l 

2(d2 -d - 1) 

d2(d+l)Vd2-l ^^^^^’-^+^ 

2 

/ I , 1 \ / T£> ^12 £d,d-|-l 


■ £12 £13 


■ £12 £34 


■ £12 £l,d-|-l 


£34 £l,d-|-l 


(d -|- l)\/d2 — 1 
d(d- 1) 

(d -|- l)\/d2 — 1 
(d- l)(3d-F2) 
d2(d + l)v/d2^ 


£l,d£l,d-|-l 


£l,d £3,d 




£d,d-|-l 


2d 


(r\J ^d,d-\-l 


^£13 £13 


£12 £34 


-3 £12 £l,d-|-l 
d 


£34 £l,d+l 


£12 £d,d+i 


£l,d £l,d-|-l 
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(2.23) 


id - 1 ) . ^ 

{d+i)Vdr^ 

id - 1) . ^ 

(d+l)^/d^ 


1 

“3 Cl es rf+i 

a 

1 

— -^^l,d^d,d+l ■ 


In the large d limit one then obtains, to leading order 

. Vd"^ — 1 


arcsin ■ 


1 


d 


+ £d,d+l + Gl,d ei,d+l + • • • 


1 


+ d + ••• - 2^M + 

+ 0(1) . 


1 2 

2 ^d,d+l 


ei2ei,d+l — ^l,d^3,d+l — ^l,d^d,d+l + ••• 


(2.24) 


To evaluate the curvature term —kJ2 ddVd -2 appearing in the gravitational lattice action one needs 
the hinge volume 14 - 2 , which is easily obtained from Eq. (2.19), by reducing d ^ d — 2, 

d^oo ~ + ••• + + ... + 612 613 + ...) + 0 ( 1 )} , (2.25) 

whereas the dehcit angle 5 is given by 

f Jrp _ 1 I 

dd = 27T- Od = 2 tt - Y S arcsin- - - h...> , (2.26) 

simplices simplices I ) 

with the expansion of the arcsin function given in Eq. (2.20). 


2.3 Evaluation of the lattice action 


We now specialize to the case where four simplices meet at a hinge. When expanded out in terms 
of the 6’s one obtains for the deficit angle 

y|- _^ 

dd = 27r — 4 .— + 'Y^ ^ ~ ^d,d+l + . . . — €i^d 6l,d+l + . . . 

simplices 


1 

d 


- 61,d — - 61 


; ^d,d+i - ^12 6i,d+i 


'-±,a 


■) + 


(2.27) 


The action contribution involving the deficit angle is then, for a single hinge. 


-k6dVd-2 


i-k) 

i-k) 


2d+2 (d_ 1) 

d\ 2+2 

2d+2 (d-l) 
dl 2+2 


1 2 



1 

jl - 2^12 + ••• 

1 {d + -- 

. ~ 6d,d+l + . . 

. ~ 6l^d 6l^d+l + . . . > 


( ~ 6d,d+l + ... — 6l^d6l^d+l + ...) . (2.28) 


It involves two types of terms: one linear in the (single) edge opposite to the hinge, as well as a 
term involving a product of two distinct edges, connecting any hinge vertex to the two vertices 



opposite to the given hinge. Since there are four simplices meeting on one hinge, one will have 4 
terms of the first type, and 4(d— 1) terms of the second type. Combining the cosmological constant 
and the curvature contributions one then obtains 


Vd 
dl 2*^/2 


^ “ 2 ^12 + ^ ^12 + 


— k ■ 2d{d — 1) ( — €d^d+i + 


~ ^l,d ^l,d+l + • • •) 

(2.29) 


The first term in the above expression refers to a single simplex, the second one to a single hinge. 
To obtain the total action, a sum over all simplices, resp. hinges, has still to be performed. 

We have also added a term ae^ in order to impose a cutoff at large edge lengths |e|. The justi¬ 
fication for this choice comes from the fact that numerical simulations show convincingly that very 
large, as well as very small, edge lengths are exponentially suppressed by the lattice gravitational 
measure, and in particular by a non-trivial interplay between the A term and the generalized trian¬ 
gle inequalities [20, 21, 24, 38, 39] (as such, a is not really a parameter that one is allowed to vary, 
and should rather be fixed to some suitable numeric value). Dropping the irrelevant constant term 
and summing over edges one obtains for the total action XJ2^d — kJ2^d Vd -2 in the large d limit 


A 




ij ^ik 


(2.30) 


up to an overall multiplicative factor y/d/d\ which will play no essential role in the following. 
The eij €ik coupling terms in the expression above can always of course be re-written in terms of 
finite differences. 


^ij (-ik — r, ^ik) T 


1 

2 "b 


■r n ^iki 


(2.31) 


and for smooth enough fields the first term on the r.h.s can be regarded as a discrete approximation 
to a derivative. 

From the action in of Eq. (2.30), one notices that its form leads naturally to a first rough 
estimate for the critical point, defined as the point where the competing A and curvature terms 
achieve comparable magnitudes, namely kc ~ \/d^. This results will be further improved below 
when we perform an explicit calculation, which takes into account the actual number of neighbors for 
each point, given a specific choice of lattice and its associated coordination number (see Eq. (2.40)). 


2.4 Action for the surface of the cross polytope 

The next step involves the choice of a specific lattice on which the action is then evaluated. One 
possibility would be the hypercubic lattice, divided into simplices as originally discussed in [18]. 

This type of lattice has 2*^—1 edges emanating from each site in d dimensions ^. Here we will evaluate 

®Which should be compared to the ~ d?/2 transverse-traceless degrees of freedom of the continuum gravitational 
field in d dimensions. The exponential growth for this particular lattice implies the existence of many redundant 
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the above action for the cross polytope fdd+i- The cross polytope I3n is the regular polytope in n 
dimensions corresponding to the convex hull of the points formed by permuting the coordinates 
(±1,0, 0,0), and has therefore 2n vertices. It is named so because its vertices are located 
equidistant from the origin, along the Cartesian axes in n-space. The cross polytope in n dimensions 
is bounded by 2^ (n — l)-simplices, has 2n vertices and 2n(n — 1) edges. In three dimensions, it 
represents the convex hull of the octahedron, while in four dimensions the cross polytope is the 
16-cell [40]. In the general case it is dual to a hypercube in n dimensions, with the ‘dual’ of a 
regular polytope being another regular polytope having one vertex in the center of each cell of the 
polytope one started with. 



Fig 1. Cross polytope /3„ with n = 8 and 2n = 16 vertices, whose surface can be used to define a simplicial 
manifold of dimension d = n — 1 = 1. For general d, the cross polytope (id+i will have 2{d -\- 1) vertices, 
connected to each other by 2d{d -\- 1) edges. 

When we consider the surface of the cross polytope in d -|- 1 dimensions, we have an object of 
dimension n—1 = d, which corresponds to a triangulated manifold with no boundary, homeomorphic 
to the sphere (as an example, see Fig. 1). The deficit angle is given to leading order by 

5d = 0 + — — {ed,d+i + 3terms -|- ei^d^i 4 +i + • • •) + 0{lld^d,e^jd) (2.32) 

degrees of freedom in the large d limit. Amusingly, it is reminiscent of the Dirac spinor case, for which the number 
of degrees of freedom is also exponential, ^ for large d. 
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and therefore close to flat in the large d limit. Indeed if the choice of triangulation is such that 
the deficit angle is not close to zero, then the discrete model leads to an average curvature whose 
magnitude is comparable to the lattice spacing or ultraviolet cutoff, which from a physical point of 
view does not seem very attractive: one obtains a spacetime with curvature radius comparable to 
the Planck length. In addition, the small fluctuation excitation spectrum for such strongly curved 
lattices looks disturbingly different from what one would expect in the continuum for transverse- 
traceless modes [21]. 

When evaluated on such a manifold the lattice action becomes 


d\ 


2 - kd^'j 1 - - 


1 /I 


1 


Dropping the 1/d correction one obtains to leading order 


dl 


2 (a - (l “ o H 4 + • • • 


(2.33) 


(2.34) 


and, up to the irrelevant constant term and an overall multiplicative factor, which can be absorbed 
into a re-scaling of the e’s, the action reduces to the simple form 






■t] 


(2.35) 


Since there are 2d{d + 1) edges in the cross polytope, one finds therefore that, at the critical point 
kd?' = A, the quadratic form in e, defined by the above action, develops 2d{d -|- 1) ~ 2d? zero 
eigenvalues 

It is worth noting here that the competing curvature {k) and cosmological constant (A) terms 
will have comparable magnitude when 

fee = ^ . (2.36) 

Here we have further allowed for the possibility that the average lattice spacing Iq = 1??)^^“^ is not 
equal to one (in other words, we have restored the appropriate overall scale for the average edge 
length, which is in fact largely determined by the value of A). This then gives for A = 1 (using the 
large-d expression for the average lattice spacing Iq, obtained later in this section in Eq. (2.38)), 
the estimate kc = \/3/(16 • 5^/^) = 0.0724 in d = 4, to be compared with kc = 0.0636(11) obtained 
in [24] by direct numerical simulation in four dimensions. Even in d = 3 one finds for A = 1, from 
Eqs. (2.36) and (2.38), kc = 2®/^/27 = 0.118, to be compared with kc = 0.112(5) obtained in [41] by 

^This result is quite close to the d^/2 zero eigenvalues expected in the continuum for large d, with the factor 
of four discrepancy presumably attributed to an underlying intrinsic ambiguity that arises when trying to identify 
lattice points with continuum points. 
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direct numerical simulation. Again, the dependence of kc on inverse powers of d is not surprising, 
as fluctuations, which are stronger in smaller dimensions, will require an increasingly larger value 
of the coupling k to make the transition happen in small dimensions. 

The average lattice spacing is easily estimated from the following argument. The volume of 
a general equilateral simplex is given by Eq. (2.7), multiplied by an additional factor of ^q. In the 
limit of small k the average volume of a simplex is largely determined by the cosmological term, 
and can therefore be computed from 


T/ ^ 1 

<y>= __log 




(2.37) 


with V{l"^) = (\/d + 1/d! 2*^/^) = CdV^- Doing the single surviving integral over P, /q°° dt^ exp(—A q/'^) 

(crfA)“^/'^r(((i + 2)/d), gives <V>= 2/dX = c^Iq- Solving this last expression for Iq then gives the 
desired result 


r — 

Lf\ — 


1 


X2/d 


2 d! 2‘^/2 


2/d 


(2.38) 


d y/d+l 

(which, for example, gives Iq = 2.153 for A = 1 in four dimensions, in reasonable agreement with the 
actual value Iq ~ 2.43 found in [24] near the transition point). The result of Eq. (2.36), extended 
to d dimensions, should then read 


kr = 






(2.39) 


which is in fact the same result as before in d = 4. Using Eq. (2.38) inserted into Eq. (2.39) one 
then obtains in the large d limit for the naturally dimensionless combination 

12/d 


kr 


T+1 


A d 


^ d^ 


m 


2 1 
d 


(2.40) 


iVd + Tl 

This result would then lead us to conclude that the above critical dimensionless ratio of couplings 
is given in the large-d limit by kc/X ~ 1/d. One should be careful though not to assign any deep 
physical significance to this result, which is only meant to help determine the critical values for the 
bare coupling constants. 

In the following we will now revert back, for simplicity, to the case of an expansion about Iq = 1. 
Returning to the partition function (and averages derived from it) associated with Eq. (2.34), we 
note that it can be formally computed via 


N 


z = 




TT 


N/2 


TT 


N/2 


i=l 


\/det M 


nili A. 


(2.41) 
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with N = 2d{d + 1). Convergence of the Gaussian integral then requires kd^ > A. From 

N 1 ^ N 1 

\nZ = — ln7 r--^lnAi ~ irlnTr- - / ds p{s) InA(s) , (2.42) 

2 ^ ^ ^ “'0 

and using the fact that for the cross polytope to leading order in l/d all eigenvalues are equal, one 
has 


logZ = 


\/d2^ 


d\ 


{k d^ — ^}j + d{d+ 1) log 


Stt/ 


^/d2i 


■+i 


d\ 


(kd? 


- A 


(2.43) 


with the first term arising from the constant term in the action, and the second term from the e-field 
Gaussian integral. Therefore the general structure, to leading order in the weak field expansion at 
large d, is log Z = ci{kd? — X) — d{d + 1) \og{kd?' — A) -|- C 2 with ci and C 2 d-dependent constants, 
and therefore d^XogZ/dX^ ~ l/{kd^ — A)^ with divergent curvature fluctuations in the vicinity of 
the critical point at kd?' = A. 


2.5 Inclusion of higher order terms 


It seems legitimate to ask what happens if the fluctuations in the e’s are large enough so that the 
quadratic approximation is no longer adequate. Then one has from Eq. (2.34), to lowest order in 
1/d, 


\/ d 2'^/2 




(2.44) 


where we have again included a cutoff term, proportional to a, for each edge. Then, again up to 
the constant term and an overall multiplicative factor, the action reduces to 

-5(A-t<^)i:4 + xT4^ (2.45) 


At strong coupling A: ^ 0, the minimum lies at a non-vanishing value of the e’s, namely eij = 
±l/y^. Since we started out with equilateral simplices with unit edges, this result is telling us 
that the edges have to be slightly extended (or shortened) to reach the minimum. As k is increased, 
the minimum eventually moves to the origin for k = A/d^. Neglecting the effects of fluctuations 
in the e fields, < e • e > — < e >^= 0, which is similar to the Landau treatment of ferromagnetic 
transitions, one then obtains 

+ ^ (2,46) 

For kd? > A the minimum is at the origin, whereas for kd? < A it moves away from it. For 
A > kd? one has a shifted minimum at cq = ±(1 — kd?/ and a total action d(eo) = — A (1 — 
kd?'IX)"^jAa. As a result cq vanishes at /c = A/d^, and so does d(eo). 
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If we apply the ideas of mean field theory, we need to keep the terms of order 1/d in Eq. (2.33). 
In the eijeik term, we assume that the fluctuations are small and replace eik by its average e. Each 
eij has 4d — 2 neighbors (edges with one vertex in common with it); this has to be divided by 2 to 
avoid double counting in the sum, so the contribution is {2d — 1) e. Then to lowest order in 1/d, 
the action is proportional to 






1 

+ 4 ^' 


E' 




(2.47) 


This gives rise to the same partition function as obtained earlier, and using it to calculate the 
average value of eij gives e, as required for consistency. 

To summarize, in this section we have developed an expansion in power of 1/d, which relies on 
a combined and simultaneous use of the weak field expansion. It can therefore be regarded as a 
double expansion in 1/d and e, valid wherever the fields are smooth enough and the geometry is 
close to flat, which presumably is the case to some extent at large distances in the vicinity of the 
lattice critical point at kc- In the next section we will develop a different and complementary 1/d 
expansion, which will not require weak fields, but will rely instead on the strong coupling (small 
k = I/SttG, or large G) limit. As such it should now be considered as a double expansion in 
1/d and k. Its validity will be in a regime where the fields are not smooth, and in fact will rely 
on considering lattice gravitational field configurations which are very far from smooth at short 
distances. 


3 Strong Coupling Expansion in Large Dimensions 


In this section we discuss the strong coupling (small k = l/(87rG)) expansion of the lattice 
gravitational partition function, first in the general case, and subsequently for large d. The resulting 
series is expected to be useful up to some k = kc, where kc is the lattice critical point (as determined 
for example from Eq. (2.40)), at which the partition function develops a singularity. It appears 
that the phase k > kc is of limited physical interest, since in that phase spacetime collapses into a 
two-dimensional manifold [20, 22, 23] (in fact, one of the first examples of compactification due to 
non-perturbative dynamics, as opposed to a specific choice of boundary conditions). 

There will be two main aspects to the following discussion. The first aspect will be the develop¬ 
ment of a systematic expansion for the partition function and the correlation functions in powers 
of k, and a number of rather general considerations that follow from it. The second main aspect 
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will be a detailed analysis and interpretation of the individual terms which appear order by order 
in the strong coupling expansion. This second part will then lead to a discussion of what happens 
for large d. 

3.1 The measure 

We will therefore first focus on the four-dimensional case, and then later exhibit its more or less 
immediate generalization to d > 4. The 4-dimensional Euclidean lattice action [16, 20, 21] contains 
the usual cosmological constant and Regge scalar curvature terms 

hatt = A ^14(^2) - kY^5h{f)Ah{f) , (3.1) 

h h 

with k = l/(87rG), and possibly additional higher derivative terms as well. The action only couples 
edges which belong either to the same simplex or to a set of neighboring simplices, and can therefore 
be considered as local, just like the continuum action. It leads to a lattice partition function defined 
as 

Ziatt = + , (3.2) 

where, as customary, the lattice ultraviolet cutoff is set equal to one (i.e. all length scales are 
measured in units of the lattice cutoff). For definiteness the measure will be of the form [20, 21, 38] 

/ idi"] = n n • (3-3) 

The lattice partition function Ziatt should be compared to the continuum Euclidean Feynman path 
integral 

Zcont = , (3.4) 

which involves a functional integration over all metrics, with functional measure [2, 42, 43] 

/ [ddi^A = /n n dgi,^{x) /n n dgi,u{x) . (3.5) 

'' X li>u X 11>U 

Since we will be doing an expansion in the kinetic term proportional to k, it will be convenient to 
include the A-term in the measure. We will set therefore in this Section 

dfi{l‘^) = [dl‘^]e-^^h^^ . (3.6) 

It should be clear that this last expression represents a fairly non-trivial quantity, both in view of 
the relative complexity of the expression for the volume of a simplex, Eq. (2.5), and because of 
the generalized triangle inequality constraints already implicit in [dl'^]. But, like the continuum 
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functional measure, it is certainly locals to the extent that each edge length appears only in the 
expression for the volume of those simplices which explicitly contain it. Also, we note that in 
general the integral / dfi can only be evaluated numerically; nevertheless this can be done, at least 
in principle, to arbitrary precision. Furthermore, A sets the overall scale and can therefore be set 
equal to one without any loss of generality (one can also conveniently normalize the integration 
measure, so that Zq = J dfj,{P) = 1, but this will not be necessary here). 

To summarize, the effective strong coupling measure of Eq. (3.6) has the properties that 1) 
it is local in the lattice metric of Eq. (2.3), to the same extent that the continuum measure is 
ultra-local, 2) it restricts all edge lengths to be positive, and 3) it imposes a soft cutoff on large 
simplices due to the A-term and the generalized triangle inequalities. Apart from these constraints, 
it does not significantly restrict the fluctuations in the lattice metric field at short distances. It will 
be the effect of the curvature term to restrict such fluctuation, by coupling the metric field between 
simplices, in the same way as the derivatives appearing in the continuum Einstein term couple the 
metric between infinitesimally close spacetime points. 

3.2 Expansion in powers of k 

Erom now on we will discuss Ziatt only, and drop the subscript latt. As a next step, Z is expanded 
in powers of k, 



It is easy to show that Z{k) = k^ is analytic at /c = 0, so this expansion is well defined, 

up to the nearest singularity in the complex k plane. An estimate for the expected location of 
such a singularity in the large-d limit was given in Eq. (2.40) of the previous section. Beyond this 
singularity Z{k) can sometimes be extended, for example, via Fade or differential approximants 
[44, 45] The above expansion is of course analogous to the high temperature expansion in 
statistical mechanics systems, where the on-site terms are treated exactly and the kinetic or hopping 
term is treated as a perturbation. Singularities in the free energy or its derivatives can usually be 
pinned down with the knowledge of a large enough number of terms in the relevant expansion 
[44]. The often surprisingly rich structure of singularities in the complex coupling plane and their 

®It is well known that a first order transition cannot affect the singularity structure of Z{k) as viewed from the 
strong coupling phase, as the free energy is Cao at a first order transition. Z{k), as defined from the strong coupling 
phase, will be non-analytic only at the second order, endpoint transition, modulo an exponentially small imaginary 
part appearing in the metastable phase, if one exists. Approaching the phase transition from the strong coupling 
side detects the physically relevant end-point singularity, where the correlation length diverges and scale invariance 
is presumably recovered [24]. 
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volume dependence has been explored in detail for some simple exactly soluble models with a finite 
number of degrees of freedom [46]. 

Next consider a fixed, arbitrary hinge on the lattice, and call the corresponding curvature term 
in the action 5A. Such a contribution will be denoted in the following, as is customary in lattice 
gauge theories, a plaquette contribution. For the average curvature on that hinge one has 

OO -1 \ ^ 

^ Ia:” [ d^i{f)5A 

< M ^^^ . (3.8) 

n=0 ■ \ h / 

After expanding out in k the resulting expression, one obtains for the cumulants 


with 


<(5A>= , 

n=0 

[ dn{f)5A 


Co = 


dn{l^ 


whereas to first order in k one has 


Cl = 


J dpi{l^)5A J dfi{l^)6A • J dfi{l^) Y^dhAh 

\ h / h 






(3.9) 


(3.10) 


(3.11) 


This last expression clearly represents a measure of the fluctuation in 5 A, namely [((X]/i 4 ^h)'^) — 
{J2h 4 Ah)‘^]/Nh, using the homogeneity properties of the lattice 5A J2h ^hAh/Nh- Equivalently, 
it can be written in an even more compact way as Nfi[{{5A)‘^) — {SA)"^]. To second order in k one 
has 


I df,{l^) M 4 Ah) J d^(/2) Y^hAh-J df^il^) M ^ 4 Ah 


C2 = 


2 / dfiil"^) 




J dfi{l^) J dfi{l^)5A J dfi{l^)5A • (j dfi{l^) Y^hAh) 


2 dfi{l^) 


+ 


dpi{V 


-<3.12) 


which now corresponds to C2 = N"^ [(((5A)^) — 3((5A)(((5A)^) + 2{5A)'^]/2. At the next order one 
has C3 = Nf^ [(((5A)^) — 4((5A)(((5A)^) — 3(((5A)^)^ + 12{{5A)^){5A)"^ — 6((5A)^]/6, and so on. Note 
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that the expressions in square parentheses become rapidly quite small, 0{1/NJ^) with increasing 
order n, as a result of large cancellations that must arise eventually between individual terms 
inside the square parentheses. In principle, a careful and systematic numerical evaluation of the 
above integrals (which is quite feasible in practice) would allow the determination of the expansion 
coefficients in k for the average curvature < 6A > to rather high order, but we shall not pursue 
this line of inquiry here 

It is advantageous to isolate in the above expressions the local fluctuation term, from those 
terms that involve correlations between different hinges. To see this, one needs to go back, for 
example, to the first order expression in Eq. (3.11) and isolate in the sum the contribution 
which contains the selected hinge with value 5A, namely 


Y^dhAh = 5A + Y.'5hAh , (3.15) 

h h 


where the primed sum indicates that the term containing SA is not included. The result is 

d^i{f){5Af (^Jdfi{f)SAy 

j Mh (^j dnil"^)^ 

E '^hAh (^J d^i{f)5^ (^J dfi{f) E 'ShAh 

One then observes the following: the first two terms describe the local fluctuation of 5A on a given 
hinge; the third and fourth terms describe correlations between 6A terms on different hinges. But 





because the action is loeal, the only non-vanishing contribution to the last two terms comes from 


edges and hinges which are in the immediate vicinity of the hinge in question. For hinges located 
®As an example, consider a non-analyticity in the average scalar curvature 


7^(fc) 


< J dx g{x) R{x) > 
< f dx ffgix) > 


(3.13) 


assumed to be of the form of an algebraic singularity at fcc, namely TZ{k) ^k^kc (kc — kY. It will lead to a 
behavior, for the general term in the series in k, of the type 


i-lffATz 


(5 - n + 1){S — n + 2)... 5 
n\ k'i-^ 


(3.14) 


Given enough terms in the series, the singularity structure can then be investigated using a variety of increasingly 
sophisticated methods [44, 47, 48, 49]. In Ref. [24] the curvature TZ{k) was computed numerically for various values of 
k, from which one can extract an approximate value for the coefficients, namely 77(fc)=—9.954 + 62.life + 195.94fe^ — 
1340.65fc® + 40483.75fc^ + 0(fc®). A better and much more accurate way would be a direct determination of the 
individual coefficients, via the edge length integrals of Eqs. (3.11) and (3.12). 
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further apart (indicated below by ^‘'notnn”) one has that their fluctuations remain uncorrelated, 
leading to a vanishing variance 


h not nn ' — 


h not nn 


dn{l^ 


dfi{l‘^ 


= 0 , 


(3.17) 


since for uncorrelated random variables X^s, < XnX^ > — < Xn >< Xm >= 0. Therefore the 
only non-vanishing contributions in the last two terms in Eq. (3.16) come from hinges which are 
close to each other. 

The above discussion makes it clear that a key quantity is the correlation between different 
plaquettes, 


<{SA)h{5A)h' > 


dix{f){5A)h{5A)h> 


(3.18) 


or, better, its connected part (denoted here by < ... >c) 


< {dA)h{6A)h'>c = < {dA)h{5A)h'> - < {6A)h >< {5A)h'> , (3.19) 


which subtracts out the trivial part of the correlation. Here again the exponentials in the numerator 
and denominator can be expanded out in powers of k, as in Eq. (3.8). The lowest order term in k 
will involve the correlation 

I d^if){5A)h{SA)h' . (3.20) 

But unless the two hinges are close to each other, they will fluctuate in an uncorrelated manner, 
with < (6 A)h ((5 A)hi > — < (5 A)^ >< {6 A)^ > = 0. In order to achieve a non-trivial correlation, 
the path between the two hinges h and h' needs to be tiled by at least as many terms from the 
product {Y.h AhT in 

I dfi{l^) {6 A)h {5 A)y 5h A^^j (3.21) 

as are needed to cover the distance I between the two hinges. One then has 

< {d A)h{6 A)h'>c ^ ~ , (3.22) 

with the correlation length ^ = \/\ log/c| ^ 0 to lowest order as A: ^ 0 (here we have used the 
usual definition of the correlation length namely that a generic correlation function is expected 

to decay as exp(—distance/^) for large separations) This last result is quite general, and holds 

^This statement, taken literally, oversimplifies the situation a bit, as depending on the spin (or tensor structure) 
of the operator appearing in the correlation function, the large distance decay of the corresponding correlator is 
determined by the lightest excitation in that specihc channel. But in the gravitational context one is mostly concerned 
with correlators involving spin two (transverse-traceless) objects. 
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for example irrespective of the boundary conditions (unless of course ^ ~ L, where L is the linear 
size of the system, in which case a path can be found which wraps around the lattice). 

But further thought reveals that the above result is in fact not completely correct, due to the 
fact that in order to achieve a non-vanishing correlation one needs, at least to lowest order, to 
connect the two hinges by a narrow tube. The previous result should then read correctly as 

< {5A)h{5A)h'>c ^ , (3.23) 

where, as will be shown in more detail below, Ud I represents (approximately) the minimal number 
of dual lattice polygons that are needed to form a closed surface connecting the hinges h and h', 
with I the actual distance (in lattice units) between the two hinges. 

3.3 Rotation matrices, Voronoi loops and closed surfaces 

Up to this point our considerations have been quite general, and therefore do not take into account 
yet the detailed nature of the local interaction expressed in the action term Ylh^hAh- It is well 
known that the deficit angle 5h describes the rotation of a vector parallel transported around 
a closed loop encircling the hinge h. This full rotation is best described in terms of a (Lorentz) 
rotation matrix R, an element of 50(4) or 50(3,1), depending on the signature of the metric, and 
whose matrix elements will depend on the specific choice of coordinates at the point in question. In 
d dimensions the corresponding objects would be SO{d) or SO{d—l, 1) rotations, in the Riemannian 
and pseudo-Riemannian case respectively 

Just as in the continuum, where the affine connection and therefore the infinitesimal rotation 
matrix is determined by the metric and its first derivatives, on the lattice the elementary rotation 
matrix between simplices is fixed by the difference between the gij's of Eq. (2.3) within 

neighboring simplices. Consider therefore a closed path T encircling a hinge h and passing through 
each of the simplices that meet at that hinge. In particular one may take T to be the boundary 
of the polyhedral dual (or Voronoi) area surrounding the hinge [20]. We recall that the Voronoi 
polyhedron dual to a vertex P is the set of all points on the lattice which are closer to P than any 
other vertex; the corresponding new vertices then represent the sites on the dual lattice. A unique 
closed parallel transport path can then be assigned to each hinge, by suitably connecting sites in 
the dual lattice. 

®The preceding observations can in fact be developed into a consistent first order (Palatini) formulation of Regge 
gravity, with suitably chosen independent transformation matrices and metrics, related to each other by a set of 
appropriate lattice equations of motion [50]. One would expect the first and second order formulations to ultimately 
describe the same quantum theory, with common universal long-distance properties. How to consistently define hnite 
rotations, frames and connections in Regge gravity was hrst discussed systematically in [51]. 
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With each neighboring pair of simplices s,s + 1 one associates a Lorentz transformation 
which describes how a given vector transforms between the local coordinate systems in these 
two simplices, 


ym 


R 


■s,s+l 




(3.24) 


The above Lorentz transformation is then directly related to the continuum path-ordered (P) 
exponential of the integral of the local affine connection T^^ via 


R 


a 

n 



path 

between simplices 
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(3.25) 


Next consider moving a vector V once around a Voronoi loop, i.e. a loop formed by Voronoi edges 
surrounding a chosen hinge. The change in the vector V is given by 


= (R - , 


(3.26) 


where R = n.R., s+i is now the total rotation matrix associated with the given hinge. Since in 
the continuum 5V is given by | > where is the antisymmetric bivector 

representing the loop area, one has the identification 


^ |3^lu^V ~ 


(R- 1) 


(3.27) 


To first order in the deficit angle 6 one then recovers the well known result 

> (3-28) 

where Ua/s represents the hinge bivector, Ua/s = 2 ^ ^af3fiu ^2 > with h and I 2 the two hinge vectors 
and Afi the area of the hinge, and use has been made of the relationship between the original 
volumes and their dual counterparts, = 2^r- As a result, one can relate the deficit angle 

directly to the effect of a complete rotation of a vector around a hinge. 


Rs,s-|-1 


u 


oShUC^') 


(3.29) 


In other words, the product of rotation matrices around the closed elementary loop describes a 
rotation in a plane perpendicular to the hinge, by an angle Sh- Equivalently, this last expression 
can be re-writtten in terms of a surface integral of the Riemann tensor, projected along the surface 
area element bivector ylp^ associated with the loop. 


[n 


R 


'5,5 + 1 


R- 


0/3 


(3.30) 
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Fig 2. Elementary closed surface tiled with parallel transport polygons. For each link of the dual lattice, the 
parallel transport matrices R are represented by an arrow. In spite of the fact that the Lorentz matrices R 
fluctuate with the local geometry, two contiguous, oppositely oriented arrows always give RR~^ = 1. 

Let us now return to the strong coupling expansion, and it will be advantageous now to focus 
on general properties of the parallel transport matrices R For smooth enough geometries, with 
small curvatures, the above rotation matrices can be chosen to be close to the identity. Small 
fluctuations in the geometry will then imply small deviations in the R’s from the identity matrix. 
But for strong coupling {k 0) it was already emphasized before that the measure / dfi{P) does 
not significantly restrict fluctuations in the lattice metric field. As a result these fields can be 
regarded in this regime as basically unconstrained random variables, only subject to the relatively 
mild constraints implicit in the measure dfi. The geometry is generally far from smooth since there 
is no coupling term to enforce long range order (the coefficient of the lattice Einstein term is zero), 
and one has as a consequence large local fluctuations in the geometry. The matrices R will therefore 
fluctuate with the local geometry, and average out to zero, or a value close to zero 

®The role of continuous rotation matrices in Regge gravity is brought out in a particularly clear way by the hrst 
order approach of Ref. [50]. 

/ cos 6 — sin 0 0 0 \ 

^°In the sense that, for example, the S'Ofd) rotation Re = ^ ^ 1 0 I out to zero when 

\ 0 0 01 / 

integrated over 9. In general an element of SO{n) is described by n(n — l)/2 independent parameters, which in 
the case at hand can be conveniently chosen as the six SO(4) Euler angles. The uniform (Haar) measure over the 
group is then d/irf(R) = d9i Jq d92 fg dOs d 64 sin 04 fjd9ssm9s dOesin^ 96- This is just a special 
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This is quite similar of course to what happens in SU{N) Yang-Mills theories, or even more 
simply in (compact) QED, where the analogs of the SO{d) rotation matrices R are phase factors 
Uf^{x) = One has there / U^{x) = 0 and / Uf^{x) Uj^{x) = 1. In addition, for two 

contiguous closed paths Ci and C 2 sharing a common side one has 

, (3.31) 

with C the slightly larger path encircling the two loops. For a closed surface tiled with many 
contiguous infinitesimal closed loops the last expression evaluates to 1, due to the divergence 
theorem. In the lattice gravity case the discrete analog of this last result is considerably more 
involved, and ultimately represents the (exact) lattice analog of the contracted Bianchi identities 
[53]. An example of a closed surface tiled with parallel transport polygons (here chosen for simplicity 
to be triangles) is shown in Fig. 2. 

We can now re-examine the question, left open earlier in this Section, of the value for the 
quantity Ud appearing in Eq. (3.23). This last quantity counts the number of polygons needed 
to obtain a closed surface around a hinge, in the framework of the strong coupling expansion for 
the curvature correlation function. For concreteness, we will consider a simplicial lattice built up 
of d-dimensional hypercubes divided up into simplices, as originally discussed in [18, 20] in the 
four-dimensional case, although similar considerations should equally apply to other semi-regular 
d-dimensional lattices as well. Simply put the issue is then; how many polygons does it take to 
form the smallest closed surface attached to two hinges, separated from each other by I lattice 
steps? 

First let us consider a slightly simpler case, namely the smallest non-trivial closed surface made 
out of elementary parallel transport loops, and built around a single given hinge. In the four¬ 
dimensional hypercubic lattice the number of triangles per edge is either 14 (for the coordinate 
edges and the hyper diagonal) or 8 (for the body and face diagonals). For a d-dimensional lattice, 
one needs the number of (d —2)-simplices on each (d —3)-simplex. This again is 14 for some (d —3)- 
simplices, and somewhat less for others. For example, using the binary notation for the vertices as 
in [18], if the vertices of the (d — 3)-simplex are taken to be (0, 0,0...), (1, 0,0,...), (1,1,0,...),... up 
to the vertex with (d — 3) Is followed by 3 Os, then the number of (d — 2)-simplices hinging on this, 
in the forward direction will be the same as the number of ways of inserting Is in the 3 remaining 
places with Os, which is 7. There will be the same number of (d — 2)-simplices in the backward 
direction. Thus for a typical (d — 3)-simplex, one needs 14 polygons to form a closed surface, 
case of the general n result [52], which reads d/irr (R) = (0"=! r(*/2)/2" 7 r"("+^)/^) n)=i with 

Q<e\< 2 'k, o<e\< TT. 
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The next step then involves considering the minimal closed surface connecting two hinges sepa¬ 
rated by I lattice steps. If one is trying to connect two polygonal half-spheres with what resembles 
a closed tube, one needs to take a path through {d — 2)-simplices connecting the (d — 3)-simplices 
at the centers of the half-spheres. Suppose the path goes through I {d — 2)-simplices, then the tube 
will consist of 26 (from the ends) plus 12(/ — 1) polygons = 12/ -|- 14. One noteworthy aspect of 
this result is that it gives a large power of k, namely ~ 12 in the notation of Eq. (3.23), but 
note that at the same time the power does not grow with d. 

In the extreme strong coupling limit this then gives, from Eqs. (3.22) and (3.23), 




k^O 


^0 

log/c^2| 


-I- . . . , 


(3.32) 


where the corrections (indicated here by the dots) arise from surfaces which are not minimal, i.e. 
deformations of the original minimal surface obtained by adding polygonal outgrowths to it, and 
therefore involving additional powers of k. 


3.4 Random surfaces and the value of the universal exponent u 


In general for k not too small the random surface spanned by the parallel transport polygons will 
have a rather complex shape. The systematic counting of such surfaces is a rather challenging task, 
say compared to a regular hypercubic lattice, in view of the simplicial nature of the underlying 
lattice geometry. When discussing the average scalar curvature, given by the expectation value of 
kd- 2 ) such a surface will be anchored on a given polyhedral loop, whereas when considering the 
correlation function of Eqs. (3.18) and (3.19) it will be anchored on two such parallel transport 
polygons, separated from each other by some fixed distance 

As one approaches the critical point, k ^ kc, one is interested in random surfaces which are of 
very large extent. Let Up be the number of polygons in the surface, and set Up = since after 
all one is describing a surface. The critical point then naturally corresponds to the appearance of 


surfaces of infinite extent. 


n„ = 


1 


kr — k 


oo 


(3.33) 


^^One might worry that the effects of large strong coupling fluctuations in the R matrices might lead to a phe¬ 
nomenon similar to confinement in non-Abelian lattice gauge theories [54, 55] . That this is most likely not the case 
can be seen from the fact that the analog of the Wilson loop W{V) (defined here as a path ordered exponential of the 
affine connection around a closed loop) does not give the static gravitational potential. The potential is instead 


determined from the correlation of (exponentials of) geodesic line segments, as in exp 




where fio is the mass of the heavy source, as discussed already in some detail in [56, 57]. The expected decay of 
near-planar Wilson loops with area A, lT(r) ~ exp(J^ R ■ ~ exp(—A/^^) ]58], where A is the minimal area 

spanned by the loop, gives instead the magnitude of the large scale, averaged curvature, operationally determined by 
the process of parallel-transporting test vectors around very large loops, and which therefore is of order R ~ 1/C^. 
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A legitimate parallel is to the simpler case of scalar field theories, where random walks of length 
T describing particle paths become of infinite extent at the critical point, situated where the 
inverse of the (renormalized) mass ^ = m~^, expressed in units of the ultraviolet cutoff, diverges 
[59, 60, 61, 62, 63], 

In the present case of polygonal random surfaces, one can provide the following concise argument 
in support of the identification in Eq. (3.33). First approximate the discrete sums over n, as they 
appear for example in the strong coupling expansion for the average curvature, Eq. (3.8) or its 
correlation, Eq. (3.18), by continuous integrals over areas 

where A = is the area of a given surface. The A'^~^ term can be regarded as counting the 
multiplicity of the surface (its entropy, in statistical mechanics terms). The exponent 7 depends 
on the specihc quantity one is looking at. For the average curvature one has 7 = —5, while for its 
derivative, the curvature fluctuation (the curvature correlation function at zero momentum), one 
expects 7 = 1 — (5. The same type of singularity is of course obtained from the original series in 
Eq. (3.34), if one assumes for the coefficients Cn 


-7 

n , 


r(l - 7 )sin 7 r(n+ 7 ) 


n ' 


TT 


1 - 


7(1-7) 

2 n 


(3.35) 


which in retrospect explains the appearance of the factor A'^~^ in Eq. (3.34). In the last step 
we have used the well-known asymptotic expansion for the binomial coefficient for large n. 
Although we know its value exactly, the integral in Eq. (3.34) can also be evaluated by standard 
saddle point methods. The saddle point is located at 

( 7 - 1 ) (7 - 1 ) kc 


A = 


log Y kc — k 


(3.36) 


Carried further, the saddle point method then leads to an approximation to the exact result for 
the quantity in Eq. (3.34), namely 

-7 


T-7 


(7 - 1 )^ ^ ^27r(7 - 1 ) (^log 


(3.37) 


which agrees with the answer given above, up to an irrelevant overall multiplicative factor. From 
this discussion one then concludes that close to the critical point very large areas dominate, as 
claimed in Eq. (3.33). 

Furthermore, one would expect that the universal geometric scaling properties of such a (closed) 
surface would not depend on its short distance details, such as whether it is constructed out 
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of say triangles or more complex polygons. In general excluded volume effects at finite d will 
provide constraints on the detailed geometry of the surface, but as d ^ oo these constraints can 
presumably be neglected and one is dealing then with a more or less unconstrained random surface. 
In the following we will assume that this is indeed the case, and that no special pathologies arise, 
such as the collapse of the random surface into narrow tube-like, lower dimensional geometric 
configurations. Then in the large d limit the problem simplifies considerably. 

Following [62], one can define the partition function for such an ensemble of unconstrained 
random surfaces as 

T r 

Zrs = [ n -/3^Aa(X) , (3.38) 

n,m=l . A 

where the integral is over d-component vectors ^n,m dehned on two-dimensional triangular lattice 
sites, with sites labeled here by integers n and m. Up to a multiplicative constant, the term 
appearing in the exponent is just the total area of the surface, written as a sum of individual 
triangle areas. Introducing the induced two-dimensional metric tensor on each triangle allows one 
to recast the above partition function in the form of a two-dimensional massless held theory, which 
in a more compact continuum notation now reads 

Zrs = const. J [dA] [d^j [dX] exp i J d^x^ {gab - Gab) ~ J , (3-39) 

with Gab = 9aX • (9;,X. The above action is now quadratic in the free massless X-helds, whose 
propagator involves A-dependent weights. We note that in the original gravitational context, the 
introduction of the coordinate vectors X(x) for describing the random surface spanned by polygons, 
originally embedded in a huctuating curved geometry, would seem plausible in view of the fact that 
as one approaches the critical point the expectation value of the scalar curvature does indeed go 
to zero [24]. 

As shown in [62], the overall size of the random surface, as embedded in the original d- 
dimensional space and suitably dehned in the discrete case as 

1 ^ 

< ^ E , (3.40) 

n,m=l 

is then immediately obtained from the free held infrared behavior of X as < X^ >~ k/T d^P/P^ ~ 
log T. Thus the mean square size of the surface increases logarithmically with the intrinsic area of 
the surface. This last result is usually interpreted as the statement that an unconstrained random 
surface has inhnite fractal (or Hausdorff) dimension. Although made of very many triangles (or 
polygons), the random surface remains quite compact in overall size, as viewed from the original 
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embedding space. In a sense, an unconstrained random surface is a much more compact object 
than an unconstrained random walk, for which < >~ T. Identifying the size of the random 

surface with the gravitational correlation length ^ then gives 


C ~ \/logr ~ I log(/Cc - . (3.41) 

k^kc 

From the definition of the exponent u, namely ^ ~ {kc — k)~^^ the above result then implies i/ = 0 
(i.e. a weak logarithmic singularity) We note that the previous result for ^ in Eq. (3.32) only 
applied to the extreme strong coupling limit A: ^ 0. 

Let us discuss next what the implications of this last result might be. As already outlined in 
Refs. [57, 64, 24], the exponent v determines the universal renormalization group evolution of the 
dimensionless coupling G = in the vicinity of the ultraviolet fixed point. In particular, 

if one defines the dimensionless function F{G) via m = = AF{G), where A is the ultraviolet 

cutoff (the inverse lattice spacing), then by differentiation of the renormalization group invariant 
quantity m, A^m(A, G'(A)) = 0, one immediately obtains the Callan-Symanzik beta function 
/3(G) [55]. From the definition 

A A (5(A) = /3(G(A)) , (3.42) 


one gets an equivalent form for the beta function in terms of the function F{G) introduced above, 
namely 


/3(G) = - 


F(G) 


(3.43) 


dF{G)/dG 

The generic beta function equation, determining the scale evolution of the coupling (obtained from 
Eq. (3.42), and identical in form to it), 


/i^G(/r) = /3(G(/r)) , (3.44) 

can then be integrated in the vicinity of the fixed point, leading to a definite relationship between the 
relevant coupling G, the renormalization group invariant (cutoff independent) quantity m = 1/^, 
and an arbitrary sliding scale fi = 1/r. Up to scales of order it determines the universal running 
of G, which will give rise to macroscopic effects provided the non-perturbative scale ^ is very large. 
In [15, 64] this scale was naturally identified with the scaled cosmological constant, which here 
would correspond to the ratio A/G. The result of Eq. (3.41) then corresponds to the limiting case 
1 / ^ 0. In the language of Refs. [15, 64], it leads in the vicinity of the fixed point to an exponentially 
small (for r/^ ^0) renormalization-group running of G(/r) or G{r), namely 

G(r) - Ge ~ e-" . (3.45) 

_ G^Gc 

four dimensions one finds for lattice quantum gravity « 1/3 instead [24, 64]. 
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All of the above was in the limit of infinite dimension. In Ref. [64] it was suggested, based on 
a simple geometric argument, that v = l/{d — \) for large d. Moreover, for the lattice theory in 
finite dimensions one finds no phase transition in d = 2 [65], u ~ 0.60 in d = 3 [41] and u ~ 0.33 in 
d = 4 [24, 64], which then leads to the (almost constant) sequence {d — 2)v = 1, 0.60 and 0.66 in the 
three cases respectively. After interpolating this last series of values with a quadratic polynomial 
in 1/d, one obtains v k, 1.9/d for large d. On the other hand in Ref. [29] the value v = l/2d was 
obtained in the same limit with a Wilson-type continuum renormalization group approach, in which 
a momentum space slicing technique is combined with a truncation to the Einstein-Hilbert action 
and a cosmological term. It seems that in either case our analytical results for the large d limit are 
consistent with, and to some extent corroborate, these previous findings. For completeness let us 
mention here that in the extreme opposite case, namely close to two dimensions, one has the by 
now well-established result v = l/{d — 2) + 0{{d — 2)^) [25, 26]. 

It is of interest to contrast the result ~ 0 for gravity in large dimensions with what one finds 
for scalar [7, 32] and gauge [66] fields, in the same limit d = oo. Known results, and what we have 
found here so far, can be combined and summarized as follows 

scalar held ^ ~ ^ 

lattice gauge held i 

lattice gravity = 0 . (3.46) 

The hrst rather well-known result is re-derived in Appendix A. The second one, obtained for non- 
Abelian gauge theories at large d, is recalled in Appendix B. It should be regarded as encouraging 
that the new value obtained here, namely 1 ^ = 0 for gravitation, appears to some extent to be 
consistent with the general trend observed for lower spin, at least at inhnite dimension. 

As far as 1/d corrections are concerned, the result obtained previously in this section hinge on 
the crucial assumption that the random surface is non-interacting, in other words that any self¬ 
intersection or folding of the surface does not carry additional statistical weights. This is similar 
to an unconstrained random walk, where the effects of path intersection and backtracking are 
neglected. While these assumptions seem legitimate at inhnite d (since there are inhnitely many 
orthogonal dimensions to move into), they are no longer valid at hnite d. As a result, the problem 
becomes much more complex, and one expects that ly will then no longer be equal to zero. Indeed 
in four dimensions z/ ~ 1/3 [24]. In the much simpler random walk case, a systematic expansion 
can be developed, leading for n intersections to an effective interaction for the scalar held 

associated with the random walk. Unfortunately in the gravitational case it is much less clear how 
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to develop such a systematic expansion. 


4 The Continuum Case 

For quantum gravity formulated in dimensions greater than four there are a number of natural 
questions that come to mind. Are there any special dimensions for gravity? How do the Feynman 
rules depend on d? What does continuum gravity look like in large dimensions? Before discussing 
the gravitational case, it might be useful to examine and contrast the somewhat simpler cases of 
scalar and vector (gauge) theories. 

4.1 Special values of d in field theories 

In scalar field theories the special role of dimension four is easily brought out by writing the action, 
simply using dimensional arguments, as 

S = \j J d^^xcpix)^ , (4.1) 

where A is the ultraviolet cutoff, Aq the bare self-coupling, mo the bare mass, and with the fields 
having canonical dimension The self-coupling is dimensionless only in dimension four, 

and above that the model is described in the long-distance, infrared limit by a free field [63]. The 
interaction term is relevant for d < 4, and irrelevant above d = 4. In particular for any d > 4 one 
can prove that the correlation length exponent n equals one half, the free field value [7, 32]. The 
long distance, infrared behavior is the same as for a free field. 

In the case of SU{N) non-Abelian gauge theories one has that the coupling is, again, dimen¬ 
sionless only in four dimensions, a well-known signature of perturbative renormalizability. Above 
four dimensions purely dimensional arguments indicate the appearance of a non-trivial ultraviolet 
hxed point (a zero of the Callan-Symanzik l3{g) function) close to the origin, 

Pia) = {d-4:)g - Pog^ + ... , (4.2) 

with a non-trivial fixed point at gl = {d — 4)//3o + 0{{d — 4)^) separating what is believed to be 
a Coulomb, non-confining phase, from the confining phase known to exist for sufficiently strong 
coupling [54]. Since the theory is not perturbatively renormalizable above four dimensions, the 
analysis of either phase is rather problematic in the continuum. The transition is characterized 
by non-trivial critical exponents, and the Green’s functions in the scaling region correspond to an 
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interacting theory, which can only be reconstructed in the Coulomb phase g < gc as an expansion 
in e = d — 4 [9]. 

One might wonder if anything special happens in dimensions d > 4, beyond what has just been 
discussed. In SU{N) Yang-Mills with (Euclidean) classical action 


one has to one loop for the divergent part of the effective action 

(1) ^ 1 26 - d g'^N 

4 - d 3 167r2 ’ 


which vanishes in d = 26, and to two loops 



1 34 (iN_y. 

4 - d 3 167r2 J ^ 


(4.3) 


(4.4) 


(4.5) 


[67]. One would be hard pressed though to conclude that the above results suggest anything 
dramatic might happen at d = 26 in the Yang-Mills case, as the change of sign in the one loop 
divergence is still counteracted by the two-(and higher-) loop terms for sufficiently large g^. It 
seems in general that the structure of the continuum theory at large d remains quite complicated 
and possibly still not amenable to a perturbative treatment. 

On the lattice on the other hand the presence of a phase transition has been clearly established 
in the large d limit, in fact largely irrespective of the specific choice of continuous symmetry group 
[66]. For the group SU{N) a critical point in g appears at 2d {2N/g'^Y = const, (with the constant 
depending of the specific choice of N), and with an exponent at the transition given hy n = 1/4 [66]. 
But it seems that finding such a transition critically hinges on using non-perturbative methods, 
which allow one to explore the strong coupling regime, and in particular the existence of two 
physically distinct phases. 

In the case of gravity, the expression analogous to Eq. (4.1) is 
A Jd'^xYgR+^ I jd^x^R^ + --- , (4.6) 


which shows the suppression of the curvature squared terms in the infrared region, by factors 
0(1/A^) when compared to the Einstein term, whose coefficient also involves a dimensionful quan¬ 
tity, namely A^“^/(167rGo) (here ao and Po, as well as Gq = A'^“^G, are taken to be dimensionless 
couplings) It then seems legitimate to ask if there are any special dimensions for gravity, in 

Adding curvature squared terms to the bare action cures the perturbative non-renormalizability problem, but 
raises new issues related to unitarity [68]. Curvature squared terms are expected to play important roles at very 
short distances, comparable to the cutoff scale, where ffuctuations in the curvature can become of order ~ A^/Go. 
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particular above d = 4. As already mentioned in the Introduction, one has d{d + l)/2 indepen¬ 
dent components of the metric in d dimensions, and the same number of algebraically independent 
components of the Ricci tensor appearing in the field equations. The contracted Bianchi identities 
reduce the count by d, and so does general coordinate invariance, leaving d(d — 3)/2 physical gravi¬ 
tational degrees of freedom in d dimensions. As a result, the number of physical degrees of freedom 
of the gravitational field grows rather rapidly (quadratically) with the number of dimensions. 

The first step is naturally to examine tree level gravity, where all loop (quantum) effects are ne¬ 
glected [1, 69, 70]. Then in the non-relativistic, static limit gravitational interactions are described 
by 


I2[T] 





□ “1 


{d - 2)-^ T/ 


d — 3 

d-2 y 


(4.7) 


where the Green’s function Q is the static limit of !/□, and = IGvrG. The above result then 
incorporates at least two well-known facts, namely that there are no Newtonian forces in d=2+l 
dimensions, and that the Einstein tensor vanishes identically in d=l+l dimensions. But nothing 
particularly noteworthy seems to happen, at least at tree level, above d = 3. At the same time, four 
spacetime dimensions is known to be the lowest dimension for which Ricci flatness does not imply 
the vanishing of the gravitational field, Rf^uXa = 0; and therefore the first dimension to allow for 
gravitational waves and their quantum counterparts, gravitons. The tree level static gravitational 
potential above d > 3 is simply obtained by Fourier transform using 




A k-x 






4yrc(/2 


-1 


(4.8) 


and therefore implies / ~ 

When quantum loop effects are turned on [3, 4], one finds that the one-loop divergence, pro¬ 
portional to curvature squared terms, vanishes on shell, 


p(i) _ 

^ div 


1 


n 


( — Ruu R^^^ + —R 
vy \ 20 ^ 120 


(4.9) 


4 — d IGtt^ 

using the well known result R^upa R^'^p^ = — R? + 4: R^y RP’^ + total derivative to eliminate 

Riemann squared terms. The complete set of one loop divergences, computed using the heat 
kernel expansion and zeta function regularization close to four dimensions, can be found in the 
comprehensive review cited in [12], and further references therein. At two loops it was shown some 
time ago [5, 6] that there is a non-removable on-shell two-loop i?^-type divergence 

1 209 G 


p(2) _ 

^ div 


4 - d 2880 (167r2)2 


J R^/^ R^ 


K,X D fd.l' 
pa 


(4.10) 
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In the last quoted reference it is argued that in the above expression the 209 arises from 11 x 19, 
with the factor of 11 coming from (26 — (i)/2, as expected from closed string theory [6]. Thus the 
latter divergence might vanish again at d = 26, but it is not expected that the same will happen 
at higher loops. 

Recent two-loop results based on the 2 -|- e expansion for gravity with a cosmological constant 
[25], inspired by the 2 -|- e of other, simpler field theory models [8, 71, 72], show the appearance of 
a non-trivial ultraviolet fixed point in the G beta function above two dimensions, 

(3{G) = (d - 2) G - ^(25 - uf) g 2 - 1^(25 - n/) + ... , (4.11) 

(for Uf massless real scalar fields minimally coupled to gravity). They could be possibly relevant 
as a first crude approximation to the four-dimensional theory (to the extent that they represent 
a manifestly gauge invariant resummation of those diagrams which can be regarded as dominant 
close to two dimensions). But unfortunately they can hardly be thought as useful in the limit 
d ^ oo, especially in view of the fact that the Borel summability in e = d — 2 [73, 74] of such an 
expansion still remains a largely open question. 

4.2 Feynman rules in d dimensions 

A direct examination of the Feynman rules for continuum gravity at large d indeed reveals the 
occurrence of some degree of simplification. But first we should clarify our conventions and notation 
for this section, which are taken from [75], and where one expands around the flat Minkowski space- 
time metric, with signature given by 77 ^,^ = diag(l, —1,—1,—1,...). The Einstein-Hilbert action in 
d dimensions is then given by 

'S'e = \/-9{x) R{x) , (4.12) 

with g{x) = and R the scalar curvature (it will also be assumed in the following that the 

bare cosmological constant is zero). Furthermore the coupling of gravity to scalar particles of mass 
m is described by the action 

5'm = ^ / d^'x \/-g{x) g^"''{x)d^(j){x)du(t){x) - (tP‘{x) . (4.13) 

Usually in perturbation theory the metric g^u{x) is expanded around the flat metric [3], by 
writing 

g^u{x) = + I^h^u{x) , (4.14) 
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with = 327rG. In the harmonic (de Bonder) gauge the graviton propagator is then given by 


, . i 'HlJipVvcr + Vpa-Viyp 

pup(j[p) — 7 : 2 i ^ 

2 + ie 


(4.15) 


which suggests that the conformal mode contribution might go away as d ^ 00 . But further 
thought reveals that this conclusion might perhaps be fallacious, as a different type of expansions 
seem to lead to slightly different conclusions. 

If one follows the method of reference [76, 77], then one defines the small fluctuation graviton 
field hfj,^{x) instead via 

{x)yj—g{x) = r]^'^ + K (x) . (4-16) 

One advantage of this expansion over the previous one is that it leads to considerably simpler 
Feynman rules, both for the graviton vertices and for the scalar-graviton vertices. A gauge fixing 
term can then be added [78, 79], for example of the form 


Kj ^ ' 


(4.17) 


as again used in [77]. The bare graviton propagator is then given simply by 

P) ( \ i VppVi^cr + gpadi'p ~ Vpiy'npcT 

D^^upaip) = X — - 

2 + ie 


(4.18) 


whose structure is now unaffected by the limit d ^ 00 . Thus with the latter definition for the 
gravitational field, there are no factors of l/{d — 2) for the graviton propagator in d dimensions; 
such factors appear instead in the expressions for the Feynman rules for the vertices. For the 
three-graviton and two ghost-graviton vertex the relevant expressions are quite complicated. The 
three-graviton vertex can be written as 


d2i 93)ai/9i,a2/32,03/93 

^(01 ^/ 3 i) {^ da 2 ( azVl 3 z ) l 32 ~ 2^“2/12^03/33) 

+9(o2'i'/32) (^^oi(o3'9/33)/3i ~ 3 X 2 ^“i/3i ^“ 3 / 33 ) 

+9(039^) (2??oi(o2^/32)/3i - 3^??oi/3ir3o2/32) 

+ ^9(o2^/32)(oi^/3i)(o39/33) + 29(q3^/93)(q2^/32)(oi9/3i) + 29(Qi^/3l)(Q3^/33)(a2 9/32) 

+9 ■ 9 ( 3 X 2 ^ 01 ( 02 ^/32)/3i 3703/33 + d^dai{a3V03)l3ida2d2 ~ ^3?oi(q23?/32)(o3 3?/33)/3i) 

+9 ■ 9 ( 3 X 23702 ( 01 3?/3i)/9237a3/33 + 3 X 23702 (03 3?/33)/32 3701 /3i ~ ‘^Va2iai'ni3i){a3V/33)l32^ 

12/^2 2 \ 

+9 '9 3 X 23703(01 3?/3i)/93?7q,2/32 + 3X237o3(o237/32)/3337oi/3i ~ 2r/Q,3(„j77,gj)(Q,237/32)/33j 

(4.19) 
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Again one notes that some terms become negligible as d ^ oo, but the remaining ones can have 
either sign, giving rise to non-trivial cancellations even for large d. The ghost-graviton vertex is 
given by 


V{ki,k2,k3)a/3,XfM 


IK 


~'>l\{aklf3)k2fi + 


(4.20) 


and the two scalar-one graviton vertex is given by 

in f 2 2 \ 

y \ Pi^iP 2 v + PiuP 2 ti - rn , 


(4.21) 


where the pi,P 2 denote the four-momenta of the incoming and outgoing scalar field, respectively. 
Finally the two scalar-two graviton vertex is given by 

2 


• 2 2 
mm 

2{d-2) 


T V/iaduX 


d-2 


d^J,l'dXa 


(4.22) 


where one pair of indices {p, v) is associated with one graviton line, and the other pair (A, a) is 
associated with the other graviton line. Again one notices some simplihcation in the limit d ^ oo. 
These rules follow readily from the expansion of the gravitational action to order ('^^)) of 
the scalar field action to order G {k^). 

The next step would involve a careful analysis of what the dominant diagrams are in the large 
d limit (still keeping in mind the serious shortcoming of assuming a vanishing bare cosmological 
constant), assuming that such a procedure remains reliable in this limit, in the sense that a com¬ 
plete resummation can be performed, and that there are no large non-perturbative, non-analytic 
contributions. But it seems so far that in the case of gravity there are conflicting claims in the 
literature [80, 81, 82] as to what exactly happens in the continuum as d ^ oo. 

In Refs. [80, 81] a gauge-invariant expansion in l/d was developed for vanishing bare cosmo¬ 
logical constant, considering both the case where the extra dimensions are non-compact and the 
case where they are highly compactihed. The observation was made that there are order-by-order 
(in 1/d) cancellations of large numbers of graphs, but the origin of such cancellation remained a 
puzzle. However, it was found that the leading term of any Green’s function was given by a set of 
disjoint bubble graphs. It was then determined that the graviton propagator acquires a physical 
pole near the Planck mass, unfortunately in a region where the validity of the expansion appears 
questionable. Finally it was claimed that at d = oo phase-space factors suppress the Feynman 
integrations and the theory is therefore finite. 

In the recent work of [82] it is also claimed that a consistent leading large-d limit exists for 
the Einstein theory without cosmological term, and that it can be constructed using a sub-class 
of planar diagrams, which seems somewhat in disagreement with the class of diagrams identified 
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in the previous references. It is then found that the large-d quantum gravity limit is well defined 
and renormalizable, provided the space-time integrations are not extended to the full d-dimensional 
space-time, in other words if the full space-time allows for compactified dimensions (the last result 
does not seem entirely surprising, as compactifying and shrinking extra dimensions leads to an 
effectively lower dimensional theory, with possibly convergent momentum integrations, depending 
on how the limits are taken). 

But it seems difficult to reconcile the above quoted results with the fact that a) the perturbative 
non-renormalizabilty issue only gets worse in the continuum as one increases the dimension, and 
b) that at least close to two dimensions an ultraviolet fixed point is known to exist, and somehow 
completely fails to show up in the large d diagrammatic treatment. The more likely scenario 
is perhaps that the theory remains perturbatively non-renormalizable even in the large d case, 
and therefore just as intractable in the continuum as the equally difficult large-d Yang-Mills case. 
So far the continuum perturbative diagrammatic treatment has not lead yet to any conclusive 
predictions about the behavior of physical gravitational observables (such as scale dependence and 
renormalization of couplings, nontrivial fixed points, anomalous scaling dimensions etc.), which 
makes it difficult to compare with recent non-perturbative lattice [24, 64] and continuum [25, 29] 
results in low dimensions, both of which project a rather different picture. 

5 Conclusions 


In this paper we have examined the lattice formulation of quantum gravity in the large d limit. 
Such a line of inquiry was stimulated by the fact that statistical systems based on local interactions 
generally tend to simplify considerably in this limit, where each point is found to be surrounded 
by a large number of neighbors, and mean field theory methods apply. Even when mean field 
theory does not apply, the hope was that the theory would simplify significantly, to a point where 
it could be solved exactly. In view of the general lack of analytical results, aside from perturbation 
theory and some other investigations restricted to low and somewhat unphysical dimensions, one 
would expect such results would help shed new light on the true non-perturbative ground state of 
quantum gravity in four dimensions. 

While d = oo at first seems rather remote from the physical case d = 4, one can make the case 
that the well known 1/N expansion of statistical mechanics system and SU{N) gauge theories (the 
planar limit) has lead to some remarkable insights into the finite N structure of these theories. 
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and in some cases even to quantitatively accurate answers for critical exponents (in the statistical 
mechanics context) and specific phenomenological predictions (for example in low energy QCD 
applications). Indeed more that once it has been argued that in the case of QCD, based on both 
theoretical and phenomenolgical arguments, that 1/N = 0 is not too remote from the physical case 
1/N = 1/3. In the same spirit, 1/d = 0 might not be as remote as it seems at hrst from the real 
world case of 1/d = 1/4 theory. 

In pursuing the 1/d expansion for gravity we have followed two somewhat complementary 
approaches. In the first approach, various terms which appear in the lattice gravitational (Regge) 
action were expanded in powers of 1/d. Since the resulting expressions are still rather cumbersome, 
we resorted to a combined weak held expansion, perturbing arbitrarily coordinated lattices built out 
of nearly equilateral simplices. The resulting expressions were then evaluated for the cross polytope, 
a triangulation of the d-dimensional sphere based on the dual of a d-dimensional hypercube. These 
were then shown to lead to a second order phase transition at a critical point kc ~ A/d, summarized 
in the result of Eq. (2.40). Near this critical point it was found that all ~ d? lattice degrees of 
freedom become massless (in the sense that all eigenvalues of the quadratic huctuation matrix 
have the same sign and approach zero), suggesting a complete disappearance of the conformal 
mode instability in the Euclidean theory at d = oo, in agreement with the naive conclusion from 
Eq. (4.15). 

The second, and perhaps more ambitious, approach was based on a combined strong coupling 
(small k = I/SttG) and large d expansion. First, in the strong coupling expansion, it was found that 
the relevant diagrams for the curvature correlation function to a given order in k can be identified 
with closed surfaces, built out of parallel transport polygons, with each polygon identified with the 
parallel transport of a test vector around an elementary loop residing within the dual lattice. We 
then argued that in the large d limit it should be possible to neglect surface self-intersections. One 
then finds that such surfaces, based on their equivalent description in terms of a two dimensional 
massless held theory, naturally give rise to a logarithmic divergence of the correlation length at the 
critical point at kc, leading in this limit to the exact (and presumably universal) result of Eq. (3.41). 

The natural question then arises, and which is difficult to ignore, of whether these large d re¬ 
sults have any relevance for a physical four-dimensional world. To the extent that the two cases 
are physically not too far apart, one would be tempted to conclude that the dependence of the 
correlation length ^ on the gravitational coupling, as expressed in Eq. (3.41), and, conversely, the 
dependence of the running gravitational coupling on as expressed in Eq. (3.45), would sug¬ 
gest, for large d, finite but exponentially small corrections to classical gravity, at least in a scaling 
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regime where the relevant distances involved are much smaller then the macroscopic curvature 
scale, r <C ^ ~ but still much larger than the Planck scale, r ^ Ip ^ ^/G. It is noteworthy 

that the quantum corrections computed here are non-analytic at r = 0, in spite of the fact that at 
short distances they become rather small, and thus provide to some extent a justification for the 
semi-classical picture of quantum gravity. In terms of the parameters relevant for vacuum struc¬ 
ture, the above-mentioned non-perturbative curvature scale then corresponds to a graviton vacuum 
condensate of order ~ 10~^^eV, extraordinarily tiny compared to the QCD color condensate 
{Aqcd = 220MeV) and the electro-weak Higgs condensate {v = 250GeH). Furthermore, as has 
been stressed before, the quantum gravity theory, at least in its present framework, does not and 
cannot provide a value for the non-perturbative curvature scale which ultimately needs to be 
fixed by phenomenological input. But, to the extent that this curvature scale clearly does not 
coincide with the Planck scale (the cutoff scale), there is some room left for it to take a very large, 
even cosmological, value. The lattice gravity model in fact provides a clear case where the naive 
identification of the curvature scale with the Planck scale can be shown to be incorrect, due to 
the highly non-trivial renormalization effects of strongly fluctuating quantum gravitational fields, 
which cleverly arrange for the two scales to differ significantly in magnitude, the more so as one 
approaches the critical point. 

Finally, in the last section, we have attempted to make contact with known results for the 
continuum theory above four dimension, and in particular those which have some degree of relevance 
for the limit d oo. Generally, and in analogy with the non-Abelian gauge theory case, it appears 
that the continuum theory does not seem to lead to the same level of simplification as the regularized 
lattice gravity model discussed in Sections 2 and 3 (and this in spite of their purported, but so 
far proven only for d = 3,4, equivalence in the lattice continuum limit). Indeed in either case 
(gravity and gauge), the issue of perturbative non-renormalizability only gets worse with increasing 
dimension. Ultimately we would tend to ascribe this state of affairs to the fact that it appears quite 
challenging to perform the needed resummation of the continuum theory with a bare cosmological 
constant (as done explicitly only close to two dimensions), perhaps an essential ingredient required 
to determine the true non-perturbative, long distance behavior of quantum gravitation - even in 
infinitely many dimensions. 
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Appendix 


A Scalar Case and Random Walks 


The scalar field case is quite straightforward and is therefore worth reproducing here. It relies 
on the well-known equivalence between the scalar field theory and the Ising model, as far as 
their critical or long distance behavior is concerned [83]. The Ising partition function is given in 
any dimension by 


Z(/3) 


^ exp ^ SiSj , 
S'i=±l _ <ij> 


(A.l) 


where <ij > denotes a sum over nearest neighbors {2d in d dimensions, for a simple cubic lattice). 
The corresponding scalar field theory is obtained by using a straightforward Gaussian integral 
representation for the Ising statistical weight, which reads 


^ ^ (2vr)-f (det/3M)-^ / d,/., e’^ , (A.2) 

Si=±l i 

and then expanding the exponent in powers of the field (j) and its derivatives. In either case the 
critical point is located where the renormalized mass of the lowest excitation vanishes. Returning 
to the Ising case, the spin susceptibility is then given by 


x(/3) = 


I 


Z{f3) 


^ ^ SoSk exp 


/3E 

<ij> 


(A.3) 


k Si=±l 

and coincides with the spin correlation function <So Sk>, summed over sites k. Equivalently, it can 
be regarded as the Fourier transform of the spin correlation function, evaluated at zero momentum. 
It is convenient to re-write the formula for the partition function as 


Z{f3) = (cosh/3)^ ^ n + , (A.4) 

5i=itl <ij> 

with t = tanh /3, N the number of sites on the lattice, and the product ranging over all links on the 
lattice. The expansion in t has an obvious diagrammatic representation [84], consisting in the case 
of x(/3) of open paths linking the site 0 to any site k, with each link appearing at most once (but 
multiple times, if the expansion in (3 is used instead). Write y = J2nXnt"', where Xn is now the 
number of open paths of length n with fixed origin. We obtain a path of length n -|- 1 by adding a 
link at its end, which can be done in 2d — 1 ways, giving Xn+i ~ ‘^dxn, and so for large d one has 
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[32] 


1 

1 - 2dt 


(A.5) 


x(/3) oc 




n=0 


Here use has been made of the fact that for large d excluded volume effects can be neglected, so 
that the factor 2d — 1 can simply be replaced by 2d. Then from x l/(p2 + m^)\p=o = (the 
spin correlation function evaluated at zero momentum) one obtains ^ ~ l/(tc — tY with tc = l/2(i 
and V = 1/2. 


B Vector case and q-coordinated Cayley Trees 


In the large d limit dominant diagrams in the strong coupling expansion of lattice gauge theories 
are represented not by surfaces, but by trees made out of three-dimensional cubes [66, 85]. In the 
case of the plaquette-plaquette correlation function, these are all the trees which can be constructed 
such that they are anchored on the two given plaquettes. 

The generating function for a q-coordinated Cayley tree [86] (a Bethe lattice with q links 
emanating from each vertex) is given by [66] 


9{t) 


u{l — I rt) 

(1 — uY 


(B.l) 


with parameters u and t related by 


t = u{l - ; (B.2) 

(7 = 3 corresponds to a trivalent or binary tree. In the SU{N) gauge case, one has t = 2d(3'^, with 
j3 ~ 2N /at strong coupling, and then the above is essentially the same as the free energy of the 
gauge theory (up to various inessential constants). Also, in the gauge case g = 6 (since a cube has 
six faces) and a new cube can be attached to any of the six faces of the original cube (again ignoring 
excluded volume effects at large d), thus creating a continuous tree made out of cubes. The free 
energy is then equal to a sum over all possible trees of arbitrary length, giving rise to hydra-like 
configurations as viewed from the diagrammatic perspective of the strong coupling expansion. 

In particular, the plaquette-plaquette correlation function is obtained from the second derivative 
of the above generating function g{t) with respect to the coupling f3. It can be represented as the 
sum of all trees of arbitrary shape (with coordination g = 6), but now with two fixed endpoints. 
Extending the analysis to general g, one can show that in fact the key result is in fact independent 
of g for g > 2. The relevant singularity in the second derivative of the free energy g{t) corresponds 
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to Uc = 1/(9 — 1). Expanding Eq. (B.2) in the vicinity of this point one finds 


t 





(B.3) 


i.e. the linear term vanishes. In the above expression tc is the critical point, 


^ (9 - 


(B.4) 


Thus tc — t ^ (u — Uc)^ for any q > 2. Eirst and second derivatives of the free energy g{t) with 
respect to t can then be calculated via 

dg ^ (9 - l)u - 1 

du {u — 1 )^ 

dg dg du 1 

dt du dt (1 — u)'? 

A ^ g(l 


which for any 9 > 2 behaves in the limit t ^ tc as 


d'^g ^ 9 (g - 2) ^ 2 “ 1 

dt'^ V2{q-1)^ Vtc-t 


(B. 6 ) 


Here use has been made of the fact that the second derivative of the free energy brings down two 
plaquette terms, giving the plaquette-plaquette correlation function, summed over both plaquette 
coordinates, and which is therefore equivalent to the Eourier transform of the plaquette-plaquette 
correlation at zero momentum. Thus one obtains the momentum space plaquette-plaquette corre¬ 
lation at zero momentum, or l/{p^ + m^)\p=o, with m = and this then gives ^ ~ l/(^c — 
and thus 1 / = 1/4 for any 9 > 2. It is further observed in [ 66 ] that the second order phase transition 
of the gauge theory described by g{t) bears a striking similarity to the condensation of branched 
polymers, with the polymer chain built out of (trees of) three-dimensional cubes. 
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